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We report the energy and centrality dependence of dynamical kurtosis 
for Au + Au collisions at ^/sWn = 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV 
at RHIC. The dynamical kurtosis of net-proton is compared to that of total- 
proton. The results are also compared with AMPT model calculations. 
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1. Introduction 

Mapping the QCD phase diagram as a function of temperature (T) and 
baryon chemical potential (/ie) is one of the main goals of current heavy 
ion experiments. The phase transition at vanishing hb is believed to be 
a crossover pQ. Model calculations suggest that the phase transition at 
large baryon chemical potential is of first order [2]. There should exist a 
point where the first-order phase transition ends, which we call the QCD 
critical point (QCP) Od]. For the infinite thermodynamical system near 
the QCP, the correlation length (£) goes to infinity and the long range 
correlations become dominant. Whereas, the correlation length can not 
be fully developed in current relativistic heavy-ion experiments due to the 
finite system size and finite evolution time. So the higher cumulant {a = 

y/((N — (N))' 2 }, S = <(Af ~Jf ))3> , k = ((iV "Jf ))4) - 3) that are more sensitive 
to correlation length are suggested [SIE]- Here, N can be the number of 
baryons in an event. Theoretical calculations of the location and property of 
the QCP meet numerical problems. In experimental aspect, the RHIC beam 
energy scan program [8] has been motivated to search for the signature of 
critical point by scan the collision energies in a wide range. It is shown that 
the fluctuations of proton number can reflect the singularity of the baryon 
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number susceptibility as expected at the QCP [9]. So the higher cumulant 
of proton number is used in current heavy ion experiments |10j . 

In experiment, the numbers of measured protons and antiprotons at 
RHIC energies are still small. In calculation of the higher cumulants, the 
statistical fluctuations are not negligible and should be subtracted from 
directly measured cumulants [TT] 

K dyn = K — K s tat- (1) 

If we have two independent Poisson distributions for proton and an- 
tiproton, respectively, then the net-proton number will follow a Skellam 
distribution, and the total-proton number will obey a new Poisson. The 
statistical contributions of kurtosis for both net-proton and total-proton 
are equal to 1/((N P ) + (Np)). They are determined only by the means of 
proton and antiproton in the selected event sample. 

By describing fluctuations of the order parameter field a near the crit- 
ical point, the calculations of the Sigma model predict that the dynamical 
kurtosis is universally negative when the critical point is approached from 
the crossover side of the phase separation line }12j . The negative dynamical 
kurtosis should be firstly observed in more peripheral collisions and sign 
change at lower incident energies. 

2. Data sample 

In this work, we use the data of Au+Au collisions at ^snn = 7.7, 11.5, 
19.6, 27, 39, 62.4 and 200 GeV collected by the STAR experiment at RHIC. 

Our analysis is restricted to particles measured within the STAR Time 
Projection Chamber (TPC) detector. The protons and antiprotons are iden- 
tified with the ionization energy loss (dE/dx) measured by the TPC within 
mid-rapidity of \y\ < 0.5 and 0.4 < pt < 0.8 GeV/c. Centrality is defined by 
the number of uncorrected charged particles within pseudorapidity \r]\ < 0.5, 
but excluding protons and antiprotons to suppress the auto-correlation ef- 
fects. For each centrality, the mean value of the number of participants is 
calculated from the Glauber model simulations. 

The centrality bin width correction is used to eliminate the volume fluc- 
tuation effect [13]. The statistical error is calculated from the Delta theo- 
rem method [2] . The systematic uncertainties are estimated by varying the 
track quality cuts and particle identification for protons and antiprotons. 

3. Results 

The energy dependence of the dynamical kurtosis of net-proton in nine 
different centrality bins is plotted in Fig. 1. The solid circles represent the 
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Fig. 1. (Color online) Dynamical kurtosis of net-proton as a function of collision 
energy in nine centrality bins. The brackets represent the systematic uncertainties. 



most peripheral (70-80%) collisions. The solid reversed triangles represent 
the most central (0-5%) collisions. We observe that below 19.6 GeV the 
dynamical kurtosis in peripheral collisions are all positive within errors. Its 
value increases towards more peripheral collisions. But above 19.6 GeV, the 
dynamical kurtosis turns to be negative in peripheral collisions. 

Figure 2 shows the centrality dependence of the dynamical kurtosis of 
net-proton in different RHIC energies. In agreement with Fig. 1, we find 
that in peripheral collisions the sign of the dynamical kurtosis of net-proton 
changes with incident energy and collision centrality. 

It is known that net-baryon number is a conserved quantity, but total- 
baryon number is not. To see whether the observed sign behavior is partic- 
ular to the conserved charge, we measure the dynamical kurtosis of total- 
proton as a function of collision energy and centrality in Fig. 3 and Fig. 4, 
respectively. The dynamical kurtosis of total-proton in peripheral collisions 
is found to be always positive at all measured energies. In central collisions, 
its value is around zero at all incident energies. Therefore, in contrary to the 
dynamical kurtosis of net-proton in peripheral collisions, we do not observe 
a sign change behavior for that of total-proton. 

In order to see if conventional non-critical effects can lead to the sign be- 
havior observed in Fig. 1, we calculate the dynamical kurtosis of net-proton 
based on the transport model of AMPT [15]. To make our calculations 
convenient for comparison with the STAR experiment, we choose the same 
kinematic cuts as used in data for protons and antiprotons. The energy 
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Fig. 2. (Color online) Centrality dependence of dynamical kurtosis of net-proton 
for Au+Au collisions from ^snn = 7.7 to 200 GeV. The shadowed areas represent 
the systematic uncertainties. 
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Fig. 3. (Color online) Dynamical kurtosis of total-proton as a function of collision 
energy in nine centrality bins. The brackets represent the systematic uncertainties. 



dependence of kurtosis from AMPT models with default version and string 
melting are shown in Fig. 5. We observe the dynamical kurtosis is positive 
in non-central collisions at all incident energies. No sign change is observed 
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Fig. 4. (Color online) Centrality dependence of dynamical kurtosis of total-proton 
for Au+Au collisions from y^sjyjv = 7.7 to 200 GeV. The shadowed areas represent 
the systematic uncertainties. 
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Fig. 5. (Color online) Dynamical kurtosis of net-proton as a function of collision 
energy from AMPT models with (a) default and (b) string melting versions. 



for two settings of the AMPT models. 



4. Summary 



The dynamical kurtosis has been measured in Au+Au collisions at y / s/v r /v 
= 7.7 to 200 GeV at STAR experiment. In peripheral collisions, the sign 
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of dynamical kurtosis of net-proton is found to change from negative to 
positive when incident energy decreases. In the contrary, the sign of that 
of total-proton in peripheral collisions keeps positive at all measured ener- 
gies. From AMPT model calculations where no critical behavior is included, 
the dynamic kurtosis for net-proton is found to be positive in non-central 
collisions for all energies. 
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